Phase diagram of spin- 1/2 bosons in one-dimensional optical lattice 
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Systems of two coupled bosonic species are studied using Mean Field Theory and Quantum Monte 
Carlo. The phase diagram is characterized both based on the mobility of the particles (Mott insu- 
lating or superfluid) and whether or not the system is magnetic (different populations for the two 
species). The phase diagram is shown to be population balanced for negative spin-dependent inter- 
actions, regardless of whether it is insulating or superfluid. For positive spin-dependent interactions, 
the superfluid phase is always polarized, the two populations are imbalanced. On the other hand, 
the Mott insulating phase with even commensurate filling has balanced populations while the odd 
commensurate filling Mott phase has balanced populations at very strong interaction and polarizes 
as the interaction gets weaker while still in the Mott phase. 
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I. INTRODUCTION 

When ultra-cold bosonic atoms with internal degrees 
of freedom, e.g. 87 Rb or 23 Na in F = 1 hyperfine state, 
are confined in magneto-optical traps, the magnetic field 
aligns the spins and the system behaves effectively as one 
of spin-0 particles. Progress in purely optical trapping 
techniques now allows the confinement of these atoms 
with their full spin degrees of freedom present 0, [|[ 
which permits the study of the interplay between mag- 
netism and superfluidity. In such systems, in addition 
to the usual spin-independent contact interaction, there 
are spin-dependent terms Q and also longer range inter- 
actions 3- Mean field calculations Q in the absence 
of the long range interactions, show that the interac- 
tions can be ferromagnetic (e.g. 87 Rb) or antifcrromag- 
netic (e.g. 23 Na), depending on the relative magnitudes 
of the scattering lengths in the singlet and quintuplet 
channels. When loaded in an optical lattice, the system 
is governed by an extended Hubbard Hamiltonian with 
spin-dependent and spin-independent terms. The phase 
diagram in the mean field and variational approximations 
has been calculated at zero and at finite tempera- 
ture 01 . 

The overall picture emerging from these calculations 
is that in the ferromagnetic case, i.e. negative spin- 
dependent interaction, the system is in a Mott insulat- 
ing (MI) phase when the filling is commensurate with 
the lattice size and the repulsive contact interaction is 
strong enough. When this interaction is weaker, or when 
the particle filling is incommensurate, the system is in a 
ferromagnetic superfluid (SF) phase. The SF-MI transi- 
tions in the ground state are all expected to be continu- 
ous. Furthermore, all Mott lobes shrink, and eventually 
disappear, as the ratio of the spin-dependent interaction 
to the spin-independent one increases. These predictions 
have been confirmed for the one-dimensional case with 
Quantum Monte Carlo (QMC) simulations 

The phase diagram for the antiferromagnctic case also 



has MI lobes at commensurate fillings and strong cou- 
pling and a SF phase at weak coupling or incommensu- 
rate fillings. However, the nature of the SF phase and 
the MI-SF transitions are different. In this case, the SF 
phase is "polar" : Superfluidity is carried either by the 
S z = component or by the S z = ±1 components. Fur- 
thermore, the Mott lobes of odd order are expected to 
be nematic in two and three dimensions and dimerized 
in one dimension ||. The MI-SF transition into these 
odd Mott lobes is generally expected to be continuous. 
In the even Mott lobes, the bosons are predicted to be 
in a singlet state which is expected to increase the sta- 
bility of the MI. This transition from non-singlet SF to 
singlet MI is argued to render the transition into the 
even Mott lobes of first order. Density Matrix Rcnor- 
malization Group (DMRG) [T^, [HI and Q MC G3 calcu " 
lations for one-dimensional systems have confirmed the 
dimerized nature of odd lobes, and determined the phase 
diagram. Further QMC work Q presented evidence of 
singlet formation in the even lobes but did not find clear 
evidence of first order phase transitions. 

A related model is obtained by coupling degenerate 
atomic ground states and degenerate excited states (see 
next section) to yield a system of two bosonic species 
sometimes referred to as spin- 1/2 bosons (HQ. As for 
the full spin- 1 case, the spin-dependent interaction can 
be ferromagnetic or antiferromagnctic. Mean field anal- 
ysis of this model [l3| yielded phase diagrams which are 
rather similar to the spin-I case. In this paper we shall 
present an extension of previous mean field treatments 
[ 1 31 ] which yields better qualitative agreement with ex- 
act QMC simulations, which we also present for the one- 
dimensional system. The paper is organized as follows: 
In Sec. II the Hamiltonian of two bosonic species is pre- 
sented. The Mean Field phase diagram is discussed in 
Sec. III. Finally, in Sec. IV Quantum Monte Carlo calcu- 
lations on one dimensional lattices verify the qualitative 
conclusions of the Mean Field theory (MFT), but provide 
quantitatively accurate values for the phase boundaries. 
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II. TWO BOSONIC SPECIES HAMILTONIAN 

We follow Ref. [l3| and consider a system of neutral 
polarizable bosonic atoms with three degenerate internal 
ground states (F g = 1) and excited states (F e = 1) char- 
acterized by the magnetic quantum number S z = 0, ±1. 
The atoms are loaded in a d-dimensional optical lattice 
produced by counterpropagating laser beams. As ex- 
plained in Ref. [l3|, in addition to generating the op- 
tical lattice, the lasers can be used to couple the internal 
ground and excited states by V and A transitions leading 
to two sets of orthogonal Bloch eigenmodes denoted re- 
spectively by and A [l3j]. The resulting Bose Hubbard 
Hamiltonian is given by 



H = -* { a l< a °3 + a U a oi) + y ^ hai ~ ^ 

+ Y C0S ^ J2 { a 0i a 0i a M a M + a M a M a 0i a 0i) 
i 

hoiUAi - IM^fiai, (1) 



where a/ ai {a ai ) creates (annihilates) a particle of "spin" 
a = 0, A on site i, Uq > is the spin- independent contact 
repulsion term and U2 is the spin-dependent interaction 
term. The number operator, h a i, counts the number of 
particles of type a on site i and /1 is the chemical po- 
tential. We take the hopping parameters t = 1 to set 
the energy scale. The phase difference, S(j) — 2(<j>o — </>a), 
gives the relative global phase between the two species. 
Equation (TIJ is the same as Eq. (9) of Ref. [l3| after one 
uses Eqs. (13) and (15) of that reference. The relative 
phase, (4>o — 4>a), was then determined by the require- 
ment that the total energy be minimized. It was argued 
in Ref. [l~3| that since U2 can be positive or negative, the 
way to minimize the energy is to have the coupling in the 
third term of Eq. (QJ be of the form — \U2\/2. This gives 
= <f>A for U-2 < and </>o = (f>\ i tt/2 for U 2 > 0. The 
Hamiltonian then becomes 



{i,j).cr cr,i 

\ u *\ / t t , t t A 

{ a 0i%i a 'Ai a Ai + a M a Ai a 0i a 0i ) 

i 

( 2 ) 

i <r,i 



However, this is not the only possible solution: The sign 
of the third term has no effect on the energy and the 
Hamiltonian obtained by putting </>o = 4>a both for U2 



positive or negative, 

H 2 = -t ( a li a *j + a tj a oi) +Y^ fl,ji ~~ ^ 

+ Y zZ { a 0i a 0i a Ai a Ai + a M a M a 0i a Oi) 
i 

+{U + J7 2 )^"oi^Ai - /i^fio-i, (3) 

i a,i 

yields exactly the same free energy as Eq. @. For the 
U2 < case, this is obvious since Eq. ^ is identical 
to Eq. J3|). For U2 > 0, this can be seen by writing 
the partition function, Z = Trc _,3ff , and noting that 
if we expand the exponential in powers of the term in 
question, only even powers will contribute. In fact, for 
U2 > 0, one can transform Eq. ([2]) into Eq. ((3j) by a simple 
phase transformation, aoi — > amc* 71 / 2 . Therefore, both 
Eqs. © and ([3]) are valid and yield the same results for 
energies and all other quantities that are invariant under 
the above phase transformation when examined with an 
exact method such as QMC. This is not necessarily so 
when approximations are made as we shall see in the 
next section. 



III. MEAN FIELD PHASE DIAGRAMS 

In this section we revisit the mean field phase diagrams 
of the model governed by Eq. (0 and . This is imple- 
mented by first writing the identity, 

al t a aj = (ch - (4)) [a aj - (O) 

+(4>^+«L(0-(4>(0- (4) 

The superfluid order parameter is {a\) = (a a ) = ip a and 
is taken to be real and uniform, i.e. independent of the 
position. The first term in Eq. Q is a small fluctuation 
and is ignored; substituting the remaining terms in Eq. [5] 
reduces H\ to the sum of single site Hamiltonians, 

H MF1 = -2id^(Vv(4 + O-0 

(J 

\U 2 \ / + t t t "\ 

2~~ V a a 0«A a A + a A a A a a ) 

+ (U + U 2 )n h A - /ti^n CT ; ( 5 ) 
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and from Eq. ([3]), H% becomes 



MF2 



-2W^(^(4 + a CT )-^ 

(T 

U 



i) 



U 2 / f t , t t 

-(U Q + U 2 )n n A - fx^hg. 



(6) 



The mean field solution is obtained by taking matrix el- 
ements of these Hamiltonians in a truncated basis of sin- 
gle site occupation number states | no 71a }, diagonaliz- 
ing the resulting matrix and then minimizing the lowest 
eigenvalue with respect to the order parameters ip a . The 
basis is truncated for large particles number (typically 
?io = riA = 12) to obtain results independent of the trun- 
cation. This gives the order parameters of the ground 
state and its eigenvector. The supcrfluid density is given 
by 



ip + ip A . 



A. Ferromagnetic: U2 < 



(7) 



The two mean field Hamiltonians, Hmfx and Hmf2 
are identical for the ferromagnetic case, U2 < 0, and their 
ground state was studied in Ref. [l3j . At commensurate 
filling and strong interaction (small t/Uo) the energy has 
one minimum at ■i/'o = V'A = corresponding to the MI. 
Figure [IJa) shows this for the Mott phase with two par- 
ticles/site, p — 2. For any non commensurate filling and 
for commensurate filling at weak coupling, there are four 
symmetric and degenerate minima with |^o| = | V-'a| 7^ 
corresponding to the supcrfluid phase, Fig. Hid). For 
later reference, we emphasize that the superfluid phase 
is symmetric in the two particle species, \ipo\ = \tpx\- 
Figures DJb,c) show how the minimum at ipo = ip\ = 
transforms continuously into four degenerate minima vis- 
ible in Fig. [Tf d) as the interaction gets weaker. The same 
behaviour is seen for the MI at all other commensurate 
fillings. 

The MF phase diagram is obtained by repeating the 
above calculation for many values of the chemical po- 
tential, p, and the interaction Uq (always keeping U2/U0 
constant). In the no-hopping limit, t/Uo — > 0, it is easy 
to see that when p satisfies the condition, 



(n - 1) (17 - ICTal) < < n (Z7 - |J7 2 |) , 



(8) 



the system is in a MI phase with n bosons per site. As 
t/Uo increases, always keeping IE/2I/E/0 constant, the MI 
region shrinks and eventually disappears giving the fa- 
miliar Mott lobes. Outside the MI the system is supcr- 
fluid. Note that U2 has the effect of shrinking the bases 
of all the Mott lobes, Fig. H When \U 2 \ = U the lobes 





-1.63 
-1.64 
-1.65 



to 



!/-'A 




FIG. 1: (Color online) The ground state energy, E, given 
by Eq. (O or ((SJ) in the ferromagnetic case, U2 = — 0.1 !7o 
at p — 2 (p/Uo = 1.3). (a) Uo — lOOi the system is in 
the MI phase, (b) Co = 25t, the system is in the MI phase 
close to the transition into the SF phase, (c) Uo = 20t, the 
system has just made the transition into the SF phase. The 
minimum has changed continuously from i/>o = ipA — to four 
degenerate minima at nonzero values of the order parameters, 
(d) Uo = 12. 5t, the system is in the superfluid phase. Note the 
degenerate symmetric minima in the SF case, \ipo\ = IV'aI 7^ 0. 
This figure agrees with Fig. 2 in Ref. [l3T |. 



disappear completely. In Fig. [2] we show mean field re- 
sults for the densities, po and p\ of the two species, the 
total density, p = pa + p\ and the corresponding super- 
fluid densities as functions of p for fixed t/Ua = 0.02 
and U2/U0 = —0.1. We see that all compressible regions, 
k = dp I dp 7^ 0, are superfluid while the incompressible 
plateaux, k = 0, are not superfluid, they are the MI. The 
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absence of discontinuous jumps in any of the densities 
as the system enters the MI indicates the transitions are 
continuous as also discussed in Ref. [l3| • We also remark 
that since in the ferromagnetic case we have |'0o| = IV'aIi 
then p S Q = p s \ and po = p\ everywhere except in the 
first Mott plateau where one of the species has zero den- 
sity, in this case pa = 0. The reason the first Mott region 
behaves this way is that in order for the two species to 
interconvert, the particles have to meet and interact. In 
the MI phase, the particles do not hop between the sites 
with the consequence that in the first MI there is only one 
particle per site and no way to interconvert the species. 
The phase diagram is shown in Fig. [3] 




U/U =-0.1, t/U =0.02, d=1 




ji/U 



o 



FIG. 2: (Color online) The densities of the two species, po 
and pa, the total density, p = p + pa, and the corresponding 
superfluid densities as functions of p/Uo. In the ferromagnetic 
case, p s o = psA and po = PA except in the first Mott lobe. 
The reason is that the particles do not interact and so there 
is no possibility to convert one species to another. 
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FIG. 3: (Color online) The phase diagram given by the mean 
field Hamiltonians, Eq. ([5]) or ([6]), in the ferromagnetic case 
U2 = — 0.1[/o. Note the shrinking of the bases of all lobes 
by \U2\/Uo. The dashed vertical line shows where the cut in 
Fig. [5] was taken. 



B. Anti-ferromagnetic: Ui > 

In the anti-ferromagnetic case, U2 > 0, the Hamiltoni- 
ans Eq. ([5]) and Eq. ^ are no longer identical. Reference 
[lH performed the mean field calculation using Eq. ([5]) 
and found that in the MI, the energy has one minimum at 
ipo = i'A = 0. As the system is taken closer to the tip of 
the Mott lobe, the energy develops local minima (Fig. 6, 
Ref. [HI); these local minima become global minima in 
the SF phase indicating that the MI-SF transition can be 
of first order. Furthermore, it was found that there is a 
continuum of degenerate minima at nonzero tpo and ij>\ 
(Fig. 6, Ref. Q). 

Doing the mean field calculation with Eq. ^ leads to 
results which arc similar in some ways to the above but 
different in very important aspects. First, for a given 
choice of t, Uq, U2 and p, we found that the two mean 
field Hamiltonians give the same ground state energy 
both in the SF and MI phases. Therefore, one cannot 
choose between them on this basis. In the MI phase, 
whether the filling is odd or even, Eq. (|6]) gives a ground 
state energy with one minimum at ipo = ipA = 0; the 
same as Eq. ([5]). As the system is brought closer to 
the lobe tip, the behavior depends on whether the filling 
in the Mott phase is odd or even. Odd filling is simi- 
lar to the ferromagnetic case: The global minimum at 
xpo = ipA = transforms continuously into four degener- 
ate minima with non-vanishing order parameter. How- 
ever, in this anti-ferromagnetic case, the minima lie on 
the tpo and ipA axes: Either -00 — while ipA 7^ 0, or 
-00 7^ while ipA = 0. This means that in the SF phase, 
the SF density is composed entirely of one species: Ei- 
ther (p s = p s0 , p s A = 0) or ( p s0 = 0, p s = p s a)- Our 
result differs from that in Ref. [13| where a continuum of 
degenerate minima was found in the SF phase. 

If the filling of the MI is even, then as the system ap- 
proaches the tip of the lobe, decreasing the interaction 
Uq, local minima develop at non zero ipo 0T fpA'- With 
Eq. ^ the local minima do not form an axisymmetric 
continuum as they do with Eq. ([5]), Ref. As for the 
odd lobes just discussed, the minima lie either on the 
ipo axis or on the ipA axis. This is illustrated in Fig. 0] 
which shows, for p = 2, the behaviour of the energy min- 
ima as the interaction is reduced taking the system from 
the MI to the SF phase. Figures E£b,c) illustrate the 
behaviour as the system is leaving the ML First, local 
minima develop with non- vanishing order parameter, in- 
dicating metastable SF, Fig. 01b). These local minima 
then become global, for lower Uq, while the previously 
global minimum at the origin becomes local indicating 
that the MI is now metastable, Fig. |4jc). Such behaviour 
signals a first order SF-MI transition for the even Mott 
phases. The nature of the transition is in agreement with 
Ref. [13| but not the nature of the SF phase. According 
to our mean field result, when the system is in the SF 
phase, superfludity is carried entirely by one species or 
the other but not by both. This result is in marked con- 
trast with that of the MFT based on Eq. §5§ and is similar 
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FIG. 5: (Color online) The densities of the two species, po 
and pa, the total density, p — po + PA, and the correspond- 
ing superfluid densities as functions of p/Uo in the anti- 
ferromagnetic case. This cut is taken just outside the third 
lobe (see Fig. [6]) so the compressibility, k = dp/dp, does not 
quite vanish. Note the discontinuous transition into the fourth 
Mott lobe. 



(c) 




FIG. 4: (Color online) The ground state energy given by 
Eq. © in the anti-ferromagnetic case, U2 = O.lUo- (a) 
Uo = 33. 3t, the global minimum is at i/jo = ipA = and 
the system is in the MI phase, (b) Uo = 12.54, the local 
minima are at non-zero tpo or tpA indicating stable MI and 
metastable SF phases, (c) Uo = 12. 2t, the global minima 
are at non-zero tpo or ipA indicating stable SF and metastable 
MI phases, (d) Uo = 11.11, four degenerate global minima 
indicating the system is in the SF phase. 

to the MFT prediction for the SF phase of the full spin-1 
model @ in the antiferromagnetic case. 

In Fig. [5] we show the densities of the two species, po 
and pA, the total density, p = p + p\, and the corre- 
sponding superfluid densities as functions of p/Uo in the 
anti-ferromagnetic case. Several features are noteworthy. 



The first, second and fourth incompressible MI plateaux 
are clearly visible but the third plateau does not quite 
form. This is because the cut shown in Fig. [5] passes 
just outside the third Mott lobe, see the phase diagram 
Fig. [6] As in the ferromagnetic case, the first Mott lobe 
is made of only one species, po m this case. In the even 
Mott lobes we have po = PA but nowhere else. In ad- 
dition, superfluidity is carried entirely by one species; in 
this case p s = p S Q and p s A = even though p\ ^ 0. 
We also note that the transition into the fourth plateau 
is discontinuous: We see clearly the discontinuous jumps 
°f Po, PA and p s as this plateau is approached. This is 
in agreement with the discussion of Fig. 2] and indicates 
that the system in Fig. Renters the fourth Mott lobe near 
its tip. This is confirmed in the full phase diagram Fig. El 

The phase diagram is obtained by calculating slices 
as in Fig. [5] for different values of p and t/Uo at fixed 
U2/U0. Figure [5] shows the mean field phase diagram for 
the case U2 = 0.1 Uq while the inset shows it for U2 = Uq. 
The boundaries of the phases do not depend much on 
whether one uses Eq. ([5]) or © for the mean field calcu- 
lation. However the nature of the SF phase depends cru- 
cially on which Hamiltonian is used as discussed above. 
For the U2/U0 = 0.1 case, and more generally for small 
U2/U0, the transition near the tip of even Mott lobes is 
first order, as discussed in connection with Figs. [5] and US 
The regions between the squares and the circles around 
the tips of the even Mott lobes in Fig. [6] denote the re- 
gions where the SF (MI) is metastable while MI (SF) is 
stable. The smaller the ratio U2/U0, the wider the first 
order transition region. On the other hand, when the 
ratio U2/U0 is large enough, there are no discontinuous 
transitions; only continuous transitions remain. The ex- 
treme case of U2 — Uq, where all odd lobes disappear, is 
shown in the inset to Fig. |B1 
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FIG. 6: (Color online) The mean field phase diagram given 
by Eq. ©for U 2 /U = 0.1. The bases of the odd MI lobes are 
reduced by U2/U0; the phase transitions near the tips of the 
even lobes are first order. The regions between the squares 
and the circles around the tips of the even lobes denote the 
regions where the SF (MI) is metastable while MI (SF) is 
stable. Inset: The phase diagram for U2/U0 = 1: all the odd 
MI lobes have vanished. 



IV. QUANTUM MONTE CARLO PHASE 
DIAGRAM 



In this section we present the results of exact Quantum 
Monte Carlo simulations for the system governed by the 
Hamiltonian Eq. @ or ([3]). As mentioned before, these 
two forms of the Hamiltonian are identical for U2 < 
and are equivalent for U2 > in that they are related by 
a simple phase rotation and do, in fact, give the same re- 
sults in QMC simulations for all the quantities we looked 
at. 

We use for our simulations the Stochastic Green Func- 
tion (SGF) algorithm [H[ with directed update This 
algorithm can be implemented both in the canonical and 
the grand canonical ensembles 17[. In this work we 



used mostly the canonical formulation where we relate 
the density to the chemical potential using 



p(N) =E(N + 1) - E(N). 



(9) 



N is the total number of particles in a system with L lat- 
tice sites and E = (H) is the total internal energy which 
is equal to the free energy in the ground state. Another 
very important quantity to characterize the phases is the 
superfluid density which is given by [l8j 



(W 2 ) 
2dtpL d - 



(10) 



in the single species case where W is the winding number 
of the bosons. In the present case, there are two species of 
particles which can be converted into each other. Conse- 
quently, the relevant winding number is the total for the 



two species and the superfluid density is then given by, 
((W + W x ) 2 ) 



2dt(3L d - 



(11) 



A. Ferromagnetic case: U2 < 

We start with ferromagnetic case, U2 < 0. As we 
did with mean field, we calculate the phase diagram in 
the (/j,/Uo,t/Uo) plane at fixed ratio U2/U0. We take 
1 / C/o = 0.1 because we want to use the same value for 
the ferro- and the anti-ferromagnetic cases and, in the 
latter case, MFT predicts first order transitions for the 
even Mott lobes for this value. With U2/U0 = —0.1, we 
study the SF-MI transition as the chemical potential (or 
density) is varied by calculating p versus p slices for dif- 
ferent fixed t/Uo- Using the canonical algorithm, this is 
done by incrementing the number of bosons by one parti- 
cle at a time, doing the simulation, and then calculating 
the chemical potential, p, using Eq. Figure [7] shows 
such a slice at t/Uo = 0.075 and clearly exhibits the first 
two incompressible Mott plateaux. Furthermore, we see 
that when the system is compressible, it is also superfluid, 
p s 7^ 0. In the Mott plateaux, the superfluid density van- 
ishes, p s = 0. 




FIG. 7: (Color online) The total particle density, p and the 
superfluid density, p a as functions of the chemical potential, 
p. The first two Mott plateaux are clearly visible. 

Next, we study the SF-MI transition at fixed commen- 
surate filling, integer p, as the coupling, t/Uo is changed. 
This exhibits the behaviour of the transition at the tips 
of the Mott lobes. In Fig. [8] we show the superfluid den- 
sity, p s versus t/Uo f° r P = li 2. We see that in both 
cases, p = 1 (Fig.|5[a)) and p = 2 (Fig.|5£b)), p s changes 
continuously as predicted by MF. 

The phase diagram is obtained by calculating p versus 
p scans for several t/Uo values is shown in Fig. [9] System 
sizes L = 20, 40, 60 were used showing very little finite 
size effects for the gap. We remark that, as expected, the 
Mott lobes shrink by Sp/Uo = 0.1 at their bases and will 



7 



l 

0.8 
,0.6 
0.4 
0.2 



(a) ■ 










□-□L=40 






^L=24 






^L=16 




I ipt w ^viptYV^Vg^ 1 P 1 1 L 








FIG. 8: (Color online) The superfluid density, p a , as a func- 
tions of t/Uo in the first (a) and second (b) Mott lobes. No 
evidence of a discontinuous jump can be seen as the system 
size increases, the transitions are continuous. 



disappear completely for |[/2| = Uq = 1. While the QMC 
and MFT phase diagrams are similar qualitatively, the 
Mott lobes given by the MFT are rounded while those 
given by exact QMC end in cusps. This is consistent 
with a critical point at the tip in the universality class of 
the two-dimensional XY model which is expected from 
general scaling arguments [2(| and first observed in QMC 
in Ref. US. 



the and A species are identical and that the peak is 
at po = p\ = 1 even though the initial state for these 
simulations was chosen to be po = 2, p\ — 0: In the 
ferromagnetic case, the density distributions, P(po) and 
P(pa), are the same both in the MI and SF phases. This 
means that the superfluid density is carried by both par- 
ticles. The same behavior was found for p = 1. This will 
be compared to the anti-ferromagnetic case below. 



P (po) = p(pa) 




t/Uo 



FIG. 10: (Color online) The density distribution of the par- 
ticles as a function of t/Uo- The total density is p = 2 and 
the values of t/Uo take the system from deep in the second 
Mott lobe to the superfluid phase. The distributions of the 
two species, and A, are identical and peak at a value of 
p/2 — 1. Compare with Fig. [15] 
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FIG. 9: (Color online) The QMC phase diagram of the ferro- 
magnetic system with U%/Uo = —0.1. All Mott lobes shrink 
by Sp/Uo = 0.1 at their bases as in the mean field case, Fig.(3] 
Note, however, that unlike the rounded lobes in the mean field 
case, the lobes here end in cusps. This is strong evidence that 
the gap vanishes exponentially, a behaviour characteristic of 
the KT universality class of the two-dimensional XY model. 
This phase diagram is very similar to that of the full spin-1 
model in Ref. Q. The dashed vertical line indicates where 
the cut in Fig. [7] was taken. 

Finally, we show in Fig. [10] the density distribution for 
total density p = 2 as t/Uo is increased taking the sys- 
tem from deep in the second Mott lobe into the superfluid 
phase as in Fig. [5fb) . We find that the distributions for 



B. Anti-ferromagnetic case: U2 > 




FIG. 11: (Color online) The total density, p as a function of 
the chemical potential, p, showing the first two incompressible 
Mott phases. Also shown is the superfluid density p s in the 
compressible phases. 

The phase diagram in the anti-ferromagnetic case, 
U2 > 0, is determined in the same way as for U% < 0: 
We calculate slices of p versus p for various values t/Uo 
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FIG. 12: (Color online) The superfluid density, p s , as a func- 
tion of t/Uo in the first (a) and second (b) Mott lobes for 
the antiferromagnetic case. No evidence of a discontinuous 
jump can be seen as the system size increases. The SF-MI 
transitions at the tips of the Mott lobes are continuous. 



keeping U2/U0 > fixed. Such a slice going through the 
first two Mott plateaux is shown in Fig. [TTJ As in the 
ferromagnetic case, we see that the compressible phase 
is superfluid, p s ^ 0. Mean field predicts that the SF- 
MI transition near the tip of the even Mott lobes is first 
order. We found no evidence for that; p s vanishes con- 
tinuously as a function of p as the MI is approached. In 
the canonical ensemble, as we have here, first order tran- 
sition are signalled by regions of negative compressibility 
[2^ . k — dp/ dp < 0. We found that k is always positive 
for all values oit/Uo indicating that the SF-MI transition 
is not first order. 




0.05 0.1 0.15 0.2 0.25 



t/u 



We confirm the nature of the SF-MI transition by 
studying the behaviour of the superfluid density as a 
function of t/Uo at fixed total density p = 1, 2. Fig- 
ure [1~2"1 shows p s versus t/Uo for the first and second Mott 
plateaux. In both cases p s changes continuously as the 
system transitions from the MI to the SF phase confirm- 
ing that all SF-MI transitions in this model arc continu- 
ous. 




FIG. 14: (Color online) The density distribution of the A 
(a) and (b) particles as functions of t/Uo. The system size 
is L = 40, f3 = 80, and the total density is p = 1. The 
values of t/Uo take the system from deep in the first Mott 
lobe to the superfluid phase. Inside the Mott lobe, the system 
is unpolarized, P(po) ~ P(pa), for t/Uo < 0.05 and both 
distributions peak at p/2 = 0.5. At t/Uo ~ 0.05, while the 
system is still in the first Mott lobe, it polarizes: The peak 
of one distribution approaches the total density, p = 1, while 
the other approaches 0. We verified that the same behaviour 
is observed for different system sizes. This transition in the 
Mott lobe is not predicted by MF. 



FIG. 13: (Color online) The phase diagram in the antiferro- 
magnetic case with data from three lattice sizes. The finite 
size effects are small. Note that, as expected, the lower bound- 
ary of the base of the second lobe shifts down by U2/U0 = 0.1 
but the upper boundary does not. When U2/U0 = 1, all odd 
order lobes disappear and only the even order ones remain. 
The vertical dashed line indicates where the p versus p cut in 
Fig. 1111 was taken. The vertical (orange) bar at t/Uo = 0.05 
in the first lobe indicates where this Mott phase polarizes. 
See the discussion of Fig. [T31 



The phase diagram is obtained, as before, by calculat- 
ing p versus p for several values of t/Uo and mapping out 
the Mott lobes. The resulting phase diagram for the first 
two lobes is shown in Fig. [T3l 

Finally, we show in Fig. Q3] that the density distribu- 
tion for total density p = 1 as t/Uo is increased taking the 
system from deep in the first Mott lobe into the super- 
fluid phase. We see that for t/Uo < 0.05, the system is 
unpolarized, the two populations have identical distribu- 
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p(pa) 




FIG. 15: (Color online) The density distribution of the A (a) 
and (b) particles as functions of t/Uo- The system size is 
L = 64, /3 — 128, and the total density is p — 2 and the values 
of t/Uo take the system from deep in the second Mott lobe 
to the superfluid phase. Inside the Mott lobe, P(po) = -P(pa) 
and both distributions peak at p/2 — 1. As soon as the 
system leaves the Mott phase and becomes superfluid, it is 
polarized: The peak of one distribution approaches the total 
density, p = 2, while the other approaches 0. The superfluid 
phase is polarized as predicted by the MFT: the superfluidity 
is carried entirely by one species. 

tions which peak at pa — pa = 0.5. At t/Uo ~ 0.05, while 
the system is still deep in the Mott lobe (see Fig. [T5|) the 
system polarizes: The population of one species drops 
abruptly to near zero while that of the other species in- 
creases to close to full filling. This polarization persists 



through the SF-MI transition. This transition inside the 
first Mott lobe is not predicted by mean field. 

In Fig. [Tn]we show the density distributions in the sec- 
ond Mott lobe. Here we find that the system polarizes 
only when it transitions from the MI to the SF phase at 
t/Uo ~ 0.17. Consequently, the SF phase in the anti- 
ferromagnetic phase is always polarized: The system is 
mostly populated by a single species and, consequently, 
superfluidity is carried by that one particle type. This 
behaviour was predicted, qualitatively, by our mean field 
calculation and is very different from the behavior in the 
ferromagnetic case. 

V. CONCLUSIONS 

In this work, we performed a mean field calculation 
of the phase diagram of the spin-1/2 bosonic Hubbard 
model, Eqs. ([2] and [3]). These models, which are derived 
from the full spin-1 model, are equivalent when treated 
exactly but we found that the mean field results depend 
on which Hamiltonian is used in the antiferromagnetic 
(U 2 > 0) case but are identical in the ferromagnetic case 
(U 2 < 0). Unlike a previous mean field result [l3| we 
found the superfluid phase is polarized when U2 > 0. 
We then used QMC to calculate the phase diagram for 
U 2 < and U 2 > 0. For U 2 < 0, we found that the 
populations are balanced both in the superfluid and in 
the incompressible phases. For U 2 > 0, we found that the 
SF phase is always polarized with one species dominating 
the population. Throughout the second Mott lobe (and 
presumably in all even Mott lobes) the populations are 
balanced. On the other hand, for t/Uo < 0.05 in the first 
Mott lobe, we found that the populations are balanced 
but that the system polarizes inside the Mott lobe at 
t/Uo ~ 0.05. This new transition inside the Mott lobe 
was not predicted before. 
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